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ABSTRACT 

This paper describes a serendipitous galaxy cluster survey that we plan to conduct with the XMM X- 
ray satellite. We have modeled the expected properties of such a survey for three different cosmological 
models, using an extended Press-Schechter (Press & Schechter 1974) formalism, combined with a detailed 
characterization of the expected capabilities of the EPIC camera on board XMM . We estimate that, over 
the ten year design lifetime of VMM, the EPIC camera will image a total of ~ 800 square degrees in 
fields suitable for the serendipitous detection of clusters of galaxies. For the presently-favored low- 
density model with a cosmological constant, our simulations predict that this survey area would yield a 
catalogue of more than 8000 clusters, ranging from poor to very rich systems, with around 750 detections 
above 2 = 1. A low-density open Universe yields similar numbers, though with a different redshift 
distribution, while a critical-density Universe gives considerably fewer clusters. This dependence of 
catalogue properties on cosmology means that the proposed survey will place strong constraints on 
the values of f^o and f^A- The survey would also facilitate a variety of follow-up projects, including 
the quantification of evolution in the cluster X-ray luminosity-temperature relation, the study of high- 
redshift galaxies via gravitational lensing, follow-up observations of the Sunyaev-Zel'dovich effect and 
foreground analyses of cosmic microwave background maps. 

Subject headings: cosmology: miscellaneous — galaxies: clusters: general — X-rays: galaxies 



1. INTRODUCTION 

Galaxy clusters are the largest gravitationally-bound 
structures in the Universe today and they are proving to 
be extremely powerful cosmological probes. In the hier- 
archical gravitational instability picture of structure for- 
mation, massive clusters arise from the extreme tail in 
the distribution of density fluctuations, so their number 
density depends critically on the cosmological parameters 
that determine the initial rms width, and the evolution 
with redshift, of that distribution. It thus follows that the 
observed cluster number density can provide strong con- 
straints on those parameters. For example, the number 
density of clusters at z = currently offers the most reli- 
able constraint (Evrard 1989; White, Efstathiou & Frenk 
1993) on the amplitude of density perturbations on small 
scales, as quantified by a& - the rms mass fluctuation in 
spheres of radius 8/i _1 Mpc, where h is the Hubble con- 
stant, Ho, in units of 100 km -1 Mpc -1 . In addition, sev- 
eral authors (White et al. 1993b; Bludman 1998; Gheller, 
Pantano & Moscardini 1998; Arnaud k Evrard 1999; Wu 
& Xue 1999) have tried to estimate the ratio of bary- 
onic to non-baryonic matter in the Universe as a whole 



from the observed baryon fraction in clusters. Perhaps the 
most exciting prospect is the possibility (Oukbir & Blan- 
chard 1992; Viana & Liddle 1996, 1999) of constraining the 
matter density, f^o, (and, perhaps, f^A = A/3JJq, where A 
is the cosmological constant) by observing the evolution 
of the number density of rich clusters. A great deal of at- 
tention (Henry 1997; Bahcall & Fan 1998; Eke et al. 1998; 
Sadat, Blanchard & Oukbir 1998; Blanchard et al. 1999; 
Borgani et al. 1999; Reichart et al. 1999a; Viana & Liddle 
1999) has been paid in recent years to this issue. To date, 
no consensus as to the value of fig has been reached, due, 
in large part, to the inadequacies of the cluster catalogues 
currently available. 

The inadequacies of current cluster catalogues moti- 
vates the creation a major new galaxy cluster catalogue 
using ESA's X-ray Multi-Mirror (XMM) satellite. The 
XMM satellite 4 was successfully launched on December 
10th 1999. It is a multi-mirror instrument, comprising 
of three Wolter type-1 X-ray telescope modules. There 
is an EPIC (European Photon Imaging Camera) imaging 
detector in the focal plane of each of the three telescope 
modules. The field of view of two of the EPIC detectors 
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is paved with 7 MOS CCDs, while the third is paved with 
12 pn CCDs. The MOS detectors share the focal plane 
of their respective telescope modules with an RGS (Re- 
flection Grating Spectrometer) camera. All five detectors 
work simultaneously, meaning that every XMM pointed 
observation will yield the type of imaging data required 
for serendipitous source detection. (This is in contrast to 
Chandra, which allows for either imaging or grating obser- 
vations, but not both at the same time.) 

To illustrate the enhanced sensitivity of XMM over other 
X-ray satellites, we have calculated, using the f akeit and 
show rates commands in XSPEC (version 10.00, Arnaud 
1996), the XMM, Chandra, ROSAT and Einstein count 
rates for an absorbed Raymond-Smith (Raymond & Smith 
1977) spectrum. The Raymond-Smith model has 4 in- 
put parameters; electron temperature (T), metallicity (Z), 
redshift (z) and normalization. For this comparison, we 
chose T = 1 keV, Z = 0.3Z Q , z = 0.1 and set the normal- 
ization so that the model spectrum had an unabsorbed flux 
of 1 x 10 -13 erg s _1 cm -2 in the 0.5-2.0 keV band. Photo- 
electric absorption (with nn — 4 x 10 20 cm -2 ) was in- 
cluded via the XSPEC wabs model, which is based on cross 
sections presented in Morrison & McCammon (1983). The 
resulting count rates are 



EPIC-pn (thin filter) 
Chandra- Acis-I 
ROSAT-PSPC 
Einstein-IPC 



0.5-10 keV 
0.5-10 keV 
0.5-2.0 keV 
0.3-3.5 keV 



0.078 s- 1 
0.017 s- 1 
0.0088 s" 
0.0040 s" 



the spatial and spectral properties of the clusters to be 
observed. In §|^ we describe how we estimated the sensi- 
tivity limits of the XCS and how we went on to use those 
limits, in combination with the results of §0, to produce 
simulated cluster catalogues. Finally, in §kT we describe 
some of the potential scientific applications of the XCS 
and discuss some of the limitations of our calculations. 

Throughout this paper we assume h = 0.5. XMM count 
rates are quoted in the 0.5-10 keV band and, except where 
stated, fluxes and luminosities are quoted in the 0.5-2.0 
keV band pass. 

2. CLUSTER SURVEYS 

Cluster catalogues have traditionally been constructed 
by identifying enhancements in the surface density of opti- 
cal galaxies on the sky (e.g Abell 1958; Abell et al. 1989) . 
While this can be made objective and algorithmic (Lums- 
den et al. 1992; Dalton et al. 1992; Postman et al. 1996), 
the projection effects that plague this approach cannot be 
overcome completely (van Haarlem et al. 1997). The small 
angular size of the X-ray emitting region in a cluster core, 
and its high contrast against the background X-ray sky, 
makes X-ray observations one of the best strategies for 
cluster detection. 

At low redshift, attention has focussed on the ROSAT 
All-Sky Survey (RASS), with a number of cluster sam- 
ples (Romer et al. 1994; Ebeling et al. 1996; Henry et 
al. 1997; Ebeling et al. 1998; Bohringer et al. 1998; De 
Grandi et al. 1999a, b) making use of its wide areal cov- 
erage. The most ambitious of the RASS surveys is the 
REFLEX survey, which covers 8235 deg 2 (com pared to 
the ~ 800 deg 2 covered by the XCS, see § |4.14j) . To date 
only a preliminary sample of REFLEX clusters has been 
published (De Grandi et al. 1999b); this sample has a flux 
limit of ~ 4 x 10~ 12 erg s _1 cm~ 2 , includes 130 clusters 
and has a maximum redshift of z = 0.308. A much larger 
sample of ~ 800 clusters, with a flux limit of 2 x 10 -12 erg 
s _1 cm~ 2 , will be released soon (Bohringer et al. 1998). 
From Table | we can see that XCS will detect a similar 
total number of clusters at z < 0.3, but this will be in a 
smaller area and to a deeper flux limit. Further, essentially 
all these clusters will be accomp anie d with serendipitous 
temperature measurements (see §5.3). By comparison, the 
largest complete sample of low-redshift cluster tempera- 
tures currently available contains only 50 objects (Blan- 
chard et al. 1999). 

At higher redshifts, the XCS will be far superior to 
the RASS-based surveys, since - with the exception of 
the NEP survey - the RASS-based surveys do not have 
the sensitivity to detect clusters beyond z ;> 0.3. The 
NEP (North Ecliptic Pole) survey has higher sensitivity 
(~ 1 x 10~ 13 erg s" 1 cm" 2 : Gioia 1998) than REFLEX, 
due to the scanning strategy of ROSAT, and has yielded 
detections of clusters as distant as z — 0.81 (Henry et 
al. 1997). Despite this enhanced sensitivity, the NEP sur- 
vey cannot compete with the XCS, since the XCS will 
cover roughly 10 times the area (800 deg 2 compared to 
84.7 deg 2 ) to roughly 10 times the depth. 

Data deeper than the RASS are, therefore, required to 
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This exercise demonstrates that XMM is ~ 4 times 
more sensitive than Chandra, ~ 10 times more sen- 
sitive than ROSAT and ~ 20 times more sensitive 
than Einstein. (The response matrices used for these 
calculations were epn_new_rmf.fits & epn_thin_arf.fits 
for XMM 5 , w215c2r_norm.rmf & w215c2r_norm.arf 
for Chandra 6 , pspcb_gain2_256.rmf for ROSAT 7 and 
ipc_90jun07_16ch.rsp for Einstein 7 .) 

The high sensitivity of XMM, combined with its wide 
field of view, excellent spatial resolution and spectral cov- 
erage, make it ideal for cluster detection out to redshifts of 
z = 1 and beyond. In this paper, we detail how an XMM 
cluster catalogue may be constructed through serendipi- 
tous detections in archival data. By examining the many 
thousands of pointing observations which will be made 
with XMM, it will be possible to build up a large sam- 
ple of clusters which extends to z ;> 1. In this paper, we 
make predictions for the numbers and types of clusters we 
hope to detect in the proposed XMM cluster survey (here- 
after XCS), and discuss the impact of the resulting cluster 
catalogue on cosmology. W e esti mate that the XCS will 
cover ~ 800 square degrees (§ 4.14 ) to an effective flux limit 
of ~ 1.5 x 10~ 14 erg s" 1 cmT z and contain (if Qq = 0.3) 



more than 8000 clusters (§ |5.3| ). 

In §|| we compare the XCS to existing and proposed 
cluster surveys. In §|| we construct a theoretical model 
for the cluster population, based on the extended Press- 
Schechter (Press & Schechter 1974) formalism of Viana & 
Liddle (1999). In §|], we describe the various assumptions 
we have made about the instrument response and about 
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detect high-redshift clusters in significant numbers, and 
several surveys have sought them through serendipitous 
detections in the fields surrounding Einstein and ROSAT 
targets. The Einstein Medium Sensitivity Survey (EMSS) 
has the largest areal coverage of any of these (734 deg 2 
above a 0.3-3.5 keV flux limit of 3.57 x fCT 12 erg s" 1 
cm~ 2 , falling to 40 deg 2 above 1.33 x 10~ 13 erg s" 1 cm~ 2 ) 
yielding a total of 104 clusters (Gioia et al. 1990; Henry 
et al. 1992). The measurement of temperatures for sev- 
eral of these clusters with redshifts above 0.3 (e.g. Henry 
1997) has led to the wide use (e.g Henry 1997; Bahcall & 
Fan 1998; Eke et al. 1998; Reichart et al. 1999a; Blanchard 
ct al. 1999; Viana & Liddle 1999; Donahue & Voit 1999) 
of the EMSS in the estimation of fio; the lack of consensus 
in the resulting constraints indicating, at least in part, the 
difficulty of using the EMSS data for such a task. 

More recently, a number of surveys (Castander et 
al. 1995; Collins et al. 1997; Jones et al 1998; Vikhlinin ct 
al. 1998a; Rosati et al. 1998; Romer et al. 2000) have been 
created from serendipitous detections in pointed ROSAT- 
PSPC observations. These surveys go much deeper than 
the EMSS, but over smaller areas: the largest single sur- 
vey is the Bright SHARC survey of Romer et al. (2000), 
which covers 179 deg 2 to a flux limit of ~ 2 x 10 -13 erg 
s _1 cm~ 2 . The deepest survey, the ROSAT Deep Cluster 
Survey, (RDCS) of Rosati et al. (1998) reaches a flux limit 
of ~ 4 x 10~ 14 erg s _1 cm~ 2 over an area of ~50 deg 2 . 

As we detail below, the XCS will detect much larger 
numbers of high-redshift clusters than the existing Ein- 
stein and ROSAT serendipitous surveys. The XCS will 
benefit not only from the increased sensitivity of XMM 
over Einstein and ROSAT (see ^), but also from XMMs 
excellent spatial and spectral resolution. These advan- 
tages would also be shared by an XMM slew survey of the 
sort proposed by Jones & Lumb (1998). However, we note 
that the relatively shallow depth of an XMM slew survey 
(~ 2 x 10~ 13 erg s _1 cm~ 2 ) means that it would detect 
few high redshift clusters, and thus have little power to 
discriminate between cosmological models. 

3. A MODEL FOR THE CLUSTER POPULATION 

Our theoretical model for the cluster population uses 
the extended Press-Schechter (Press & Schechter 1974) 
formalism of Viana & Liddle (1999). The validity of this 
general approach has been demonstrated by comparison 
with iV-body simulations (Eke, Cole & Frenk 1996; Col- 
berg et al. 1998; Tormen 1998). We refer the reader to 
Viana & Liddle (1999) for a more detailed description of 
the method. 

We consider herein three cosmological models, namely: 

1. the currently-favored spatially-flat, low-density cos- 
mology with Sl = 0.3 and ^a = 0.7, 

2. the Einstein - de Sitter critical density cosmology, 
with fio = 1 (and Qa = 0) 

3. an open cosmology, with fio = 0.3 (Oa = 0). 

We note that throughout this paper, we calculate lumi- 
nosity and angular diameter distances as follows. For the 
two Ha = models, we use the standard exact form due to 
Mattig (Mattig 1958). For the model with non-zero Q\, 
we use the approximate form derived by Pen (1999). This 



form is perfectly adequate for our purposes, since it has 
an error of <J 1% for flat cosmologies. 

We assume that structure formation proceeds through 
gravitational instability from a Gaussian distribution of 
primordial density perturbations with a scale-invariant 
power spectrum. The extended Press-Schechter formalism 
enables us to compute the number density of clusters as a 
function of redshift. The version we use (Viana & Liddle 
1999) includes a tracking of the merger histories of clus- 
ters in order to account properly for their time of formation 
when relating their mass to their temperature. The mass 
to temperature conversion is normalized so as to reproduce 
the results from the hydrodynamical simulations of White 
et al. 1993b and Bryan & Norman 1998, with the exten- 
sion to open or flat cosmologies with an arbitrary value for 
Oq performed using the expressions given in Viana & Lid- 
dle (1996). In the following we have only included systems 
with T > 2 keV, because the Press-Schechter formalism 
becomes unreliable at low temperatures. 

There is a weak dependence of Press-Schechter results 
on the current shape of the linear power spectrum of den- 
sity fluctuations, so, for definiteness, we have used a Cold 
Dark Matter (CDM) power spectrum with shape param- 
eter T (Efstathiou, Bond & White 1992) equal to 0.23, as 
suggested by some analyses of galaxy clustering (Peacock 
& Dodds 1994; Viana & Liddle 1996; but see Mann, Pea- 
cock & Heavens 1998). The normalization of the power 
spectrum is that of Viana & Liddle (1999), ensuring that 
the present-day abundance of high-temperature clusters is 
recovered. Similar cluster-based normalizations were also 
obtained by Eke, Cole & Frenk 1996, Pen 1998, Borgani et 
al. 1999, Blanchard et al. 1999 and Henry 2000. Further, 
these models give a good fit to the COBE 4-year data 
(e.g. Tegmark 1996). 

In Fig. [l| we plot the cumulative number, N , of clus- 
ters with temperatures greater than 2, 4, and 6 keV in 
the whole sky as function of flux cut, /, for the three cos- 
mologies. To derive the N(f) functions from the Press- 
Schechter results we have to assume a conversion from 
cluster temperature to luminosity. For this we use the 
empirical cluster luminosity-temperature relation derived 
by Allen & Fabian (1998, AF98); 

T = 1.66i£ 429 (1) 

where the temperature, T, is in keV and the bolometric 
luminosity, Lx, in is units of 10 ergs . Observations 
to date present no evidence for significant evolution of 
the L x -T relation out to z ~ 0.4 (AF98; Mushotzky 
& Scharf 1997; Reichart, Castander & Nichol 1999b), but 
nothing is known beyond that. For the purposes of our 
calculations, we assume that equation (|l|) holds at all red- 
shifts, but stress (as discussed further in §|^) that one of 
the principal scientific results of the XCS will be a greatly 
improved understanding of the Lx — T relation and its 
evolution with redshift. 

In Press-Schechter theory, the relative abundance of 
galaxy clusters of a given mass at two given redshifts de- 
pends only on the growth rate of perturbations, which in 
turn depends only on fi and Fig. [l] shows that N(f) 
varies between the three cosmologies more dramatically 
as temperature increases (note the different scales in the 
three panels). Moreover, below T ~ 4 keV, it is possible 
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that the mass-temperature relation has been significantly 
influenced by heat injection into the intergalactic medium. 
For these reasons, we will largely focus our discussion on 
clusters with X-ray temperatures in excess of 4 keV (or lu- 
minosities £ 2.6 x 10 44 erg s" 1 , based on the AF98 L x -T 
relation). Although we note that, in practice, the opti- 
mum (i.e. the one that minimizes the errors on cosmolog- 
ical parameter estimates) temperature limit for the XCS 
will probably not be exactly 4 keV. 

From Fig. [j] it is clear that distinguishing between high 
and low values of flo is relatively straight forward, but 
that discriminating between open and flat models with 
the same value of fio = 0.3 is much harder. To do so, one 
needs to have access to clusters at sufficiently high redshift, 
as illustrated by Fig. || which shows the cumulative flux 
distribution of clusters at z > 1: the predicted numbers of 
T > 6 keV clusters at z > 1 in the two Slo = 0.3 models 
differ by as much as a factor of four at faint flux limits, 
compared to less than than a factor of two for z > 0. In 
Fig. [3] we plot the analogous curves for clusters at z < 0.3, 
to demonstrate that cluster catalogues limited to low red- 
shift (such as those produced by the various RASS-based 
surveys or that to be produced by the Sloan Digital Sky 
Survey of Gunn et al. 1998) are very poor at constraining 
cosmological parameters: it is necessary to reach z ;> 0.5 
to get a sufficiently long lever arm in cosmological time for 
the sensitivity of the growth rate of density perturbations 
to cosmology to become apparent. 

The curves plotted in Fig. || are flat for fluxes fainter 
than ~ 10~ 13 erg s _1 cm~ 2 , indicating that all z < 0.3 
clusters with temperatures above T = 2 keV can be de- 
tected above that flux limit. This is no surprise, of course: 
by assumption (through equation Q) clusters above a cer- 
tain temperature also exceed a certain luminosity, so there 
clearly must be a flux level at which they are all visible 
if a redshift limit is imposed. What is more interesting is 
that this asymptoting behavior is also seen in the curves in 
Fig. [j], for which there is no redshift limit. The levelling off 
of the curves in Fig. pi results from the fact that clusters do 
not exist at arbitrarily high redshifts in a hierarchical uni- 
verse: it takes a certain amount of time to accumulate the 
matter making up a cluster of a given mass (and, hence, 
temperature and luminosity). So, if one can reach a suffi- 
ciently faint flux limit, one can look along one's past light- 
cone beyond the epoch when the first cluster of a particular 
mass was formed. For the high temperature (T > 4 keV) 
clusters important for cosmological parameter estimation, 
this asymptote is being reached at a depth (~ 10~ 14 erg 
s~ 4 cm~ 2 ) wh ich is comparable to that to be reached by 
the XCS (§5.3). This implies that there would be no point 
in ever performing a deeper survey than the XCS if the 
sole purpose of that survey was to detect T > 4 keV clus- 
ters; only a survey with a wider sky coverage would yield 
more detections (and hence tighter constraints on f2o and 
Oa)- Although we note that, by going deeper, one can 
obtain useful information about individual clusters, such 
as their spatial morphology and their temperature profile. 

4. ASSUMPTIONS AND SIMPLIFICATIONS MADE WHEN 
SIMULATING THE SENSITIVITY LIMITS OF THE XCS 

In we predicted the number of clusters on the sky, 
N, as a function of X-ray flux, /. If clusters were the 
only sources in the X-ray sky, then N(f) would describe 



the cluster catalogue resulting from the idealized situation 
of an all-sky survey performed by an instrument with no 
internal background, a vanishingly narrow point spread 
function (so that confusion noise is zero) and a uniform 
flux limit. In reality, the XCS will have a non-uniform 
flux limit, will only cover a fraction of the sky and will 
have to detect clusters against a significant X-ray back- 
ground. Moreover, since the X-ray sky is dominated by 
point sources (e.g. AGN), a crucial step in its construc- 
tion will be the differentiation of point-like from extended 
sources. 

So, in order to predict how many, and what type of, 
clusters will be included in the XCS catalog, we need to 
make various assumptions about the sensitivity and spec- 
tral response of the instrument; about the surface bright- 
ness profiles and spectral properties of the clusters we ex- 
pect to observe; and about the properties of the X-ray 
background. We detail these, and other, assumptions be- 
low. 

4.1. EPIC-pn only 

We have simplified our calculations by concentrating 
only on the EPIC-pn camera. This is because the EPIC- 
MOS cameras receive only 50% of the flux from their re- 
spective telescope modules (the other 50% in each goes 
to an RGS), and because the MOS CCDs are intrinsically 
less sensitive than the pn CCDs. An additional advantage 
is that the simulation of the catalogue selection function 
will also be simplified if data from only one camera are 
used. However, once the clusters have been detected, the 
EPIC-MOS data can be used to help parameterize the 
cluster morphology and spectrum. By prudent use of the 
EPIC-MOS data, the percentage of clusters with accom- 
panying temperatures should increase over that suggested 
by Table |. 

4.2. Minimum Detection Threshold of 8a 

We adopt a minimum detection threshold of 8a. This 
is because the XCS will have to rely on source extent to 
differentiate clusters from point-like X-ray sources, such as 
stars and AGN, and it has been shown (e.g. Greg Wirth, 
private communication) that extent measures can only be 
derived with confidence for sources detected at ;> 8a. The 
Bright SHARC survey (Romer et al. 2000) also used a 
minimum detection threshold of 8a for this reason. 

4.3. Detection Significance Computed using Inner 50% of 

Flux 

When calculating the detection significance, we only 
consider the inner 50% of the total cluster flux. (We define 
the radius of the region enclosing this flux as r§Q.) This 
is because automated source detection algorithms tend to 
underestimate the count rates of extended sources. For 
example, the wavelet transform method adopted by the 
Bright SHARC survey, to analyse ROSAT data, underes- 
timates the count rate of z > 0.15 clusters by a factor of 
2.1 (Romer et al. 2000). This is a conservative assump- 
tion: we would hope that more efficient cluster selection 
algorithms will be developed, to make full use of the higher 
quality XMM data. 
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4.4. Clusters Follow Spherically-Symmetric Isothermal 
(3=2/3 Model 

To estimate r$o values, we assume that all clusters can 
be modeled as spherically-symmetric systems that follow 
an isothermal /3-profile 



A Raymond-Smith spectrum with an unabsorbed flux of 
1 x 10 



I = 



1 + {r/r^YP- 1 / 2 



(2) 



where I is the surface brightness at radius r, r c is the 
core radius and —3/3 is the asymptotic radial fall-off of 
the ICM (intracluster medium) density distribution. We 
adopt (3 = 2/3 throughout, since this is a typical value for 
rich clusters (Jones & Forman 1984, Mohr et al. 1999), 
although we note that its value for any given cluster can 
vary in the range 0.4 <, f3 <, 0.9. For a cluster described 

by equation (^), (3 = 2/3 gives rso = V3r c . We adopt 
this model for the cluster surface brightness, as it has been 
shown (e.g. Mohr et al. 1999) to describes the azimuthally- 
averaged cluster emission in the ROSAT band- pass (0.5- 
2.0 keV) very well. However, as we discuss in § |6.2.3j , the 
use of such a simplistic model is one of the major limita- 
tions of our calculations. 

4.5. Non-Evolving Core Radius - Luminosity Relation 

We further assume that the core radius follows the rela- 
tion 



0.2 



250 /L 44 . , 
h^{— ] kpC ' 



(3) 



where L44 is the rest-frame luminosity in the 0.5-2.0 keV 
band in units of 10 44 erg s -1 . This relation was proposed 
by Jones et al. (1998) and has been shown to agree with 
measured values of r c for clusters with luminosities in the 
range 10 43 erg s _1 to 10 45 erg s _1 . We assume that the 
core radius does not evolve (as shown by Vikhlinin et 
al. 1998b). A better understanding of the luminosity - 
core radius relation (in particular, whether it evolves with 
redshift) should result from forthcoming observations of 
known clusters with Chandra and XMM . 

4.6. Cluster Count Rates 

To determine how the XMM count rate varies with clus- 
ter parameters, we used the fake it and show rates com- 
mands in XSPEC and assumed that the cluster X-ray emis- 
sion can be described by absorbed Raymond-Smith spec- 
tra. We calculated unabsorbed fluxes and on-axis XMM 
count rates for spectra with 12 different temperatures 
(1 < T < 12 keV, in 1 keV increments) and 40 different 
redshifts (0.05 < z < 2.0, in Az = 0.05 increments); 480 
spectra in all. Th rough out we kept the metallicity fixed at 
Z = 0.3Z Q (see §4.10), the normalization fixed at 1 and 
the Galacti c hi column density fixed at n# = 4 x 10 20 



cm (see §4.11). The ignore command was used to limit 



the count rate calculation to the 0.5-10 keV band. (The 
full energy range over which EPIC-pn is sensitive is 0.1-11 
keV.) 

These 480 calculations provided us with the count rate 
to flux conversion factors that were used to define the sur- 
vey sensitivity limits in §|^. To illustrate how these conver- 
sion factors vary with T and z, we provide some examples: 

8 astro . estec . esa.nl/XMM/user/uhb_top .html 
9 Available from legacy.gsfc.nasa.gov 



~ 13 erg s -1 cm -2 will yield 0.078 EPIC-pn counts 
s" 1 (0.5-10 keV) when T = 1 keV and z = 0.1. A spec- 
trum with the same flux will yield 0.071 counts s _1 when 
T = 1 keV and z = 1, 0.100 counts s" 1 when T = 10 keV 
and z = 0.1, and 0.093 counts s _1 when T = 10 keV and 
z=l. 

We note that the AF98 Lx — T relation used in sec- 
tion §^ was constructed using a slightly different plasma 
model to that used here; mekal (Kaastra & Mewe 1993) 
in XSPEC rather than raymond. However, this does not 
present a problem for this study since we limit our discus- 
sion to T > 2 keV clusters (the predictions of mekal and 
raymond are very similar above T J> 1 keV). 

4.7. Cosmic Background Count Rate 

We calculated the cosmic background using a model that 
includes two thermal Galactic components (modeled with 
absorbed Raymond-Smith spectra) and a power-law ex- 
tragalactic component. The first thermal component had 
a temperature of 0.0258 keV, a metallicity of Zq, a red- 
shift of z = 0, a normalization of 2.5 x 10~ 6 and a hy- 
drogen column density of nn = 1 x 10 17 cm -2 (Labov & 
Bowycr 1991). The second thermal component had a tem- 
perature of 0.0947 keV, a metallicity of Zq, a redshift of 
z = 0, a normalization of 3.0 x 10~ 6 and a hydrogen col- 
umn density of nn = 6x 10 19 cm~ 2 (Rocchia et al. 1984). 
The power-law component had an index of a = 1.4, a 
normalization of 9.32 x 10~ 7 (Chen et al. 1997) and a 
hydrogen column density of nn — 4 x 10 20 cm -2 . The 
adopted cosmic background model yields a count rate of 
2.6 x 10~ 3 s _1 arcmin -2 in the 0.5-10 keV band. We note 
that, when calculating the signal-to-noise of cluster detec- 
tions, we adjust the cosmic backgr oun d count rate by the 
appropriate vignetting factor (see §4.9). The true external 
background (made up of solar, Galactic and extragalactic 
components) is known to vary considerably across the sky, 
but most of this variation is confined to low energies (< 1 
keV, Snowden et al. 1997), and so this should not have a 
significant effect on the average signal-to-noise values for 
cluster detections we calculate in the 0.5-10 keV band. 



4.8. Particle Background Count Rate 

We calculated the particle background using the ex- 
pected internal background rates quoted in the XMM 
Users' Handbook?: 3.0 x 10 -4 counts cm -2 s -1 keV -1 for 
the EPIC-pn detector. The spectrum of the internal back- 
ground is expected to be flat, so the integrated count rate 
in the 0.5-10 keV band is 9.5 x 3.0 x 10 -4 counts cm -2 
s . We then converted from cm -2 to arcmin -2 to obtain 
a rate of 1.4 x 10 -4 counts s -1 arcmin -2 (4.1" corresponds 
to 150 /mi at the detector). 

4.9. Vignetting Correction 
The count rate to flux conversion factors calculated us- 



ing XSPEC ( §|4.6| ) refer to the on-axis response of the EPIC- 
pn. In order to account for how these conversion factors 
vary with off-axis angle, we had to calculate vignetting 
corrections. We did this as follows; using the QUICKSIM 9 
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package written by Steve Snowden, we created fake EPIC- 
pn images of a point source, with a Raymond-Smith spec- 
trum, in the absence of particle and cosmic backgrounds. 
By placing the source at various places in the field of 
view, we were able to measure how the count rate var- 
ied as a function of off-axis angle. The vignetting fac- 
tor changes smoothly across the field of view, so we de- 
cided to break up the field of view into five 3' wide annuli 
(9 = 1'.5,4'.5,7'.5,10'.5,13'.5 respectively). For each an- 
nulus we calculated the mean vignetting factor for a point 
source with a T = 4 keV spectrum; this was found to be 
0.987, 0.892, 0.734, 0.578 and 0.520 respectively. We used 
a single temperature for this calculation because we found 
the vignetting factor to be essentially independent of tem- 
perature; the on-axis sensitivity is 2.09 times that of the 
sensitivity at 9 — 12' for a T — 1 keV spectrum, compared 
to 2.14 for a T — 8 keV spectrum. We also made a mega- 
second QUICKSIM simulation of the cosmic background, in 
the absence of sources and a particle background, to con- 
firm that these average vignetting factors also apply to the 
cosmic background. 

4.f0. Constant ICM Metallicity: Z = O.3Z 

The X-ray emission from an astrophysical plasma is a 
function of its metallicity. For example, we calculate that 
for a Raymond-Smith spectrum with an unabsorbed flux 
of 1 x 10 -13 erg s _1 cm" 2 , the count rate varies from 0.077 
s- 1 to 0.078 s- 1 to 0.079 s -1 for Z = Q.lZp>, 0.3Z o and 



Z Q respectively (T = 1 keV, z = 0.1, n H = 4x 10 20 cm -2 ). 
Because of this weak dependence of count rate to metal- 
licity, we adopt a constant value of Z = O.3Z , as this 
is typical of rich clusters: Fukazawa et al. (1998) found 
that the ensemble-averaged iron abundance was 0.3 ±0.02 
based on ASCA observations of 40 nearby clusters of galax- 
ies. Further, we assume that metallicity does not evolve; 
up to z ~ 1 there is observational support for this from 
Tsuru et al. (1997) and Schindler (1999), and theoretical 
support from calculations by Martinelli et al. (2000). 

4.11. Constant hi Column Density: nn = 4 x 10 20 cm~ 2 

Neutral hydrogen gas along the line of sight towards 
a cluster, particularly within our own galaxy, absorbs a 
large fraction of the emitted X-rays at low (<; 0.5 keV) 
energies. Since we do not know what the actual distribu- 
tion of hydrogen column densities will be in the XCS, we 
have adopted a single value, nn = 4 x 10 20 cm -2 , which 
is typical for high Galactic latitudes. 

The effect of column density on count rates is not large. 
For Raymond-Smith spectra with unabsorbed fluxes of 
1 x 10~ 13 erg s -1 cm~ 2 , the count rate varies from 0.086 
s^ 1 to 0.078 s^ 1 to 0.064 s" 1 for n H equal to 1, 4 and 10 
x 10 20 cm' 2 respectively (Z = 0.3Z Q , T = 1 keV, z = 0.1). 
Our adoption of 4 x 10 20 cm -2 is on the conservative side; 
many of the regions explored by XMM will have lower 
nn values. For example, of the 37 clusters in the Bright 
SHARC survey (Romer et al. 2000), all but 9 were detected 
in regions with n# < 4 x 10 20 cm~ 2 . This means that the 
number of clusters eventually detected by the XCS could 
well be higher than suggested by Table |j. 

4.12. EPIC Thin Filter used for all Observations 



The EPIC cameras are sensitive not only to X-rays, 
but also to optical photons. Optical blocking filters (thin, 
medium or thick) are used to minimize the number of pho- 
tons entering the detector. For our calculations we use the 
response functions corresponding to the thin filter only. 

The choice of optical filter has an even smaller effect 
on count rates than column density. For a Raymond- 
Smith spectrum with an unabsorbed flux of 1 x 10 -13 



erg s 1 cm 2 , 



the count rates vary from 0.078 s 



to 



0.076 s _1 to 0.06 s' 1 when the thin, medium and thick 
filters are respectively in place [Z = 0.3Z Q , T = 1 keV, 
z = 0.1, nn = 4 x 10 20 cm -2 ). To calculate the count rate 
through the thin, medium and thick filters we used the files 
epn_thin_arf.fits, pn_med_arf.fits and epn_thick_arf.fits 10 
respectively. It is unlikely that any XMM pointings that 
require the thick filter (i.e. those with bright stars in their 
field of view) will be suitable for serendipitous cluster de- 
tection and so we can safely discount the effects of filter 
choice on the cluster numbers presented in Table ||. 

4.13. Bolometric and K-corrections 

We calculated K-corrections and bolometric corrections 
using xspec. The K-correction was defined as the ratio 
of the unabsorbed flux in the observed energy band to the 
unabsorbed flux in the redshifted energy band 



Kln-hi — 



r 



hi(l+z) 
(1+*) 



f v du 



(4) 



where lo and hi are the limits of the observed energy band, 
e.g. 0.5 and 10.0 keV. When calculating K-corrections, 
the redshift of each Raymond-Smith spectrum was set to 
2 = 0. Quadratic fits to the K-corrections, as a function 
of (1 + z) were derived for each of the input temperatures 
(1 to 12 keV); 



Kio-hi = c + 6(1 + z) + o(l + zf 



(5) 



where a, 6, c are the coefficients of the fits, see Table 0. 

The bolometric correction was defined as the ratio of the 
unabsorbed flux in a pseudo-bolometric band of 0.01-50 
keV to the unabsorbed flux in the observed energy band, 
i.e.; 



Bin-] 



Ii h o f» dv 



(6) 



Setting the redshift of the Raymond-Smith spectrum to 
z = 0, B values were calculated for each of the 12 input 
temperatures. The bolometric corrections are listed in Ta- 
ble for the 0.5-2.0 keV and 0.5-10 keV energy bands. 

To illustrate how the bolometric and K-corrections were 
applied, we provide an example. Consider a cluster with 
temperature T = 4 keV, a redshift of z — 1, and an 
unabsorbed flux in the 0.5-2.0 keV (observed) band of 
1 x 10~ 13 ergs -1 cm~ 2 . The K-corrected flux in the 0.5- 
2.0 keV (rest frame) band is 0.836 x 10~ 13 ergs -1 cm -2 
(from equation || and Table |l|) . The bolometric flux for 
this cluster is then 3.04 times this, from Table |[ 



5 All available from astro.estec.esa.nl in the directory tree /pub/XMM/EPIC/March99/RESP0NSES 
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4.14. Exposure Time Distribution and Area of the Survey 

In £|5.3 l w e combine our sensitivity limit calculations 
(§ |5.2| fc §p.l|) with our model cluster population (§EjJ) in 
order to predict the properties of the XCS. To do so re- 
quires us to assume both an areal coverage and an expo- 
sure time distribution for the survey. We do not know 
what the exposure time distribution will be for the thou- 
sands of pointings that will eventually comprise the XMM 
archive. So, for the purposes of this paper, we assume that 
the exposure times will be distributed in the same way as 
they are for 760 pointings in the XMM Guaranteed Time 
Observations (GTO, Table |). The 760 GTO pointings 
have exposure times that range from 5 ks to 95 ks, with 
an average of 22.3 ks. For comparison, we also list in Ta- 
ble H the distribution of the exposure times in the XMM 
A01 program. 

For the areal coverage, we use a total value of 800 deg 2 
(as justified below). How ever, we note that our treatment 
of vignetting effects (§ [4.9| ) forces us to break this total area 
up into five bins when creating mock cluster catalogues. 
These bins correspond to the five adopted off- axis annuli, 
which cover 4.3%, 13.1%, 21.7%, 30.4% and 30.2% of the 
total area respectively. 

The EPIC field of view covers a 30' diameter circle and 
the CCD arrangement of the pn camera provides an ac- 
tive area of 649 square arcminutes. If EPIC operates for 
the full ten years of the XMM design lifetime, and XMM 
makes an average of three pointings per day, then the total 
area imaged by EPIC will be ~ 2000 deg 2 . (Although the 
average exposure time of the 760 GTO pointings is 22.3 
kilo-seconds, or 3.9 pointings per day, overheads, such as 
the ~ 5 ks telescope settling time, mean that three point- 
ings per day is a more realistic estimate.) Unfortunately, 
not all of the ~ 2000 deg 2 will be available for building 
serendipitous cluster catalogues. Experience from ROSAT 
suggests that only ~ 40% of pointings are likely to be suit- 
able, the rest being either at low Galactic latitude, over- 
lapping previously studied fields, or have pointing targets 
extending over most of the field of view. Therefore, we 
estimate that the XCS will cover ~ 800 deg 2 . 

5. SIMULATION RESULTS 

The mechanisms outlined in §|] allow us to simulate 
the sensi tivit y of the XCS in terms of bot h c luster de- 
tection (§5.1) and temperature estimation (§5.2). The re- 
sults of these simulations, when combined with the Press- 
Schechter predictions descri bed in §[|, allow us to predict 
the properties of the XCS (§[5.3|). 

5.1. Calculation of Sensitivity Limits for Detection 



As stated in §4.2, our criterion for source detection is 



that it should be made with count statistics significant 
at the 8cr level at least, so that it is possible to deter- 
mine whether the source is extended or not. For the XCS 
predictions, we have calculated the bolometric luminosity 
that would yield an 8a detection (I/so-) for each of 144,000 
different parameter combinations. These 144,000 combi- 
nations comprise of 3 cosmologies (§||), 12 temperatures 
(1 < T < 12 keV, in 1 keV increments), 40 redshifts 
(0.05 < z < 2, in Az = 0.05 increments), 5 off-axis an- 
gles (0 = 1'.5,4.'5,7'.5,10'.5,13'.5, [Q and 20 exposure 
times (5 < t < 100 ks, in 5 ks increments). 



We determine the 144,000 L% a values iteratively as fol- 
lows. For a particular Oa, T, z, and t comb inat ion, 
we start by calculating the half-flux radius, r$o (§ [4.3| ), in 
arcminutes for a given input r c value (r$o — VSr c , §fOj). 
Next, we c alcu late the total number of background counts, 
that would fall in a circle of radius 7-50 in 



N, (M §4. 



the exposure time t (where N takes into account the ef- 
fects of vignetting on the cosmic background at off-axis 
angle 9). Once N is known, we can calculate the number 
of cluster counts, S, that would need to fall inside r$o to 
yield S/N > 8. By multiplying S by 2/t, we obtain the 
corresponding total count rate (where the factor of 2 ac- 
counts for those photons lying outside the r§o radius and 
the factor of 1/t converts from total counts to a count 
rate). This count rate can then be converted into a flux 



from equation 
luminosities dif 



by dividing by the appropriate vignetting factor (§4.9) and 
then multiplyin g by the appropriate count rate to flux con- 
version, factor (§4.6). Using the appropriate K-correction 
(§ [4.13 ), we calculate the corresponding (0.5-2.0 keV) lu- 
minosity. We compare this luminosity to the one obtained 
using the input value of r c . If the two 
er by more than 20%, we recalculate r c 
and repeat the whole procedure. We always start the it- 
eration with r c = 250 kpc. Usually the process converges 
after only 1 or 2 adjustments to r c . After convergence, 
w e defi ne L% a using the appropriate bolometric correction 

(§P§- 

We note that, across all 144,000 calculations, the small- 
est rso value used to derive an L 8cr value was ~ 20". To 
confirm that clusters of this size, and larger, will be flagged 
as extended sources, we have computed, using QUICKSIM, 
how the radius enclosing half the flux in a model XMM 
point spread function varies with off-axis angle. We find 
that the maximum size of this radius is only ~ 13". There- 
fore, it follows that all clusters detected at > 8a will be 
flagged as extended sources, i.e. the completeness of the 
XCS should not suffer by the imposition of an extent cri- 
terion. 

It would be impractical to provide tables listing the re- 
sults of all 144,000 calculations, so instead we give a few 
illustrative examples. These examples are based on calcu- 
lations in the first radial bin (9 — 1'.5) and in an fio = 1, 
Oa = cosmology. 

1. The highest number of counts required inside r§Q for 
an 8 a detection was 255, arising for a cluster at red- 
shift z = 2 with a temperature T = 1 keV, observed 
for 100 ks. Such a cluster would have a total count 
rate of 0.0051 count s" 1 , a flux of 7.06 xlO -15 erg 
s _1 cm~ 2 and a luminosity of 7.48 xlO 44 erg s _1 . 

2. The lowest number of counts required inside r$Q for 
an 8 a detection was 77, arising for a cluster at red- 
shift z = 0.35 with a temperature of T = 12 keV, 
observed for 5 ks. Such a cluster would have a total 
count rate of 0.0308 count s^ 1 , a flux of 3.11 xlO" 14 
erg s~ x cm~ 2 and a luminosity of 0.16 xlO 44 erg s _1 . 

3. All clusters lying at z < 1.5 which are brighter than 

will be detected at > 8a in a 5 ks pointing. At 
z = 1.5 an L+ cluster has a flux of 4.24 xl0~ 14 erg 
s _1 cm~ 2 . Here we assume that L*, the 'knee' in the 
Schechter function fit to the X-ray cluster luminos- 
ity function, has a value of 4.8 x 10 44 erg s _1 : this 
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is the average of the values found by De Grandi et 
al. (1999a) and Ebeling et al. (1998). Such a cluster 
has a temperature of 5.5 keV based on the Lx — T 
relation of AF98, equation (fy). 



4. In the average exposure time (22.3 ks, § 4.14 ) of the 
760 GTO pointings, it will be possible to detect clus- 
ters brighter than 2.4 x 10 44 erg s _1 (i.e. 0.5 x L+) 
out to redshifts of z ~ 1.6. (Such a cluster has a 
temperature of 3.7 keV based on the Lx — T relation 
of AF98.) 

5.2. Calculation of Sensitivity Limits for Temperature 
Estimation 

The EPIC-pn camera is able to estimate the energies 
of all incident photons, so it can perform low-resolution 
spectroscopy as well as broad-band imaging. This means 
that we are able to estimate temperatures for the clus- 
ters we detect, provided the signal-to-noise ratio of the 
source spectrum is sufficiently high. This is clearly a great 
advantage, since, for those clusters for which it is possi- 
ble, we shall not need to obtain follow-up observations 
to determine their temperature (c.f. the ASCA follow-up 
of EMSS clusters by Henry 1997). Cluster temperature 
measurements are important for cosmological parameter 
estimation, since T is more readily related to cluster mass 
(in terms of which theoretical predictions are made) than 
is (the more easily measured) X-ray luminosity. 

To assess the extent to which we will be able to take 
advantage of the XMM spectral resolution, we have calcu- 
lated the minimum bolometric luminosity (Lp) that would 
yield a temperature estimate for each of our 144,000 pa- 
rameter combinations. We stress that, in most cases, de- 
termination of the redshift of the cluster will be required 
before its temperature can be estimated. This remains 
true even for spectra of high signal-to-noise, due to the de- 
generacy between temperature and redshift in the spectral 
fitting when thermal bremsstrahlung is the dominant emis- 
sion process. However, if there is significant line emission 
in addition to the bremsstrahlung radiation (which is espe- 
cially true for low temperature, high metallicity plasmas), 
it is sometimes possible to measure the redshifts from emis- 
sion features such as the 7 keV Fe Line. Mushotzky (1994) 
and others have shown that this technique works and we 
will certainly apply it where possible to XCS data. Alter- 
natively, we might also expect to be able to obtain crude 
redshift estimates using the measured flux and extent: this 
method is cosmology-dependent, but is still worthy of fur- 
ther investigation. 

Given that the redshift will be known prior to the spec- 
tral fitting, it will be possible to choose the metric aperture 
size most suitable for temperature measurements for each 
cluster. To reflect this, we allow the radius of apertures 
used in our temperature sensitivity limit calculations to 
vary (with the constraint that it must never be smaller 
than rso), so as to include the maximum number of pho- 
tons but without being swamped by the background. (This 
is in contrast to our detection sensitivity calculations, for 
which we always used r§o , since the detection software will 
most likely only pick out the central ~ 50% of the cluster 
flux, see §h3). We also set the additional criteria that the 
number of background counts in the aperture must never 
exceed the number of cluster counts, and that the cluster 



counts must never be less than 1000. The aperture sizes 
thus chosen varied from to rgg and the number of back- 
ground counts inside these apertures varied from ~250 to 
~1500 (with an average value of ~600). 

The accuracy to which temperatures can be estimated 
depends on three factors: the cluster redshift, the cluster 
temperature, and the signal-to-noise of the spectrum. We 
illustrate this via Figure ^, which shows the input versus 
fitted temperatures for representative values of the clus- 
ter redshift, cluster temperature and background count 
rate. For this Figure, we created, and then fitted, 20 fake 
spectra (with nn = 4 x 10 20 cm~ 2 and Z = 0.3Z Q ) for 
each of the listed temperature-redshift-background com- 
binations using the XSPEC commands fakeit and fit re- 
spectively. The mean and standard deviation of the 20 
fits are plotted in Figure [|. These fits were all performed 
on spectra containing 1000 counts, because only about 1% 
of the 144,000 calculations produced background counts 
- and hence cluster counts - that exceeded 1000. From 
the Figure, it is clear that a spectrum of 1000 counts will 
yield temperatures estimates of varying accuracies, with 
the most accurate values being derived for the lowest tem- 
perature systems. The lowest accuracy results will come 
from high-temperature clusters at low redshift with high 
background count rates. We note that the accuracy im- 
proves with redshift for high-temperature systems because 
the 'knee' in the thermal bremstrahllung spectrum moves 
to lower energies, where XMM has more effective area. 

Some examples of our temperature sensitivity limits are 
as follows (assuming 8 = 1'.5, do = 1 and fl\ = 0). Cool 
clusters (T — 2 keV) will yield temperature measurements 
only out to z ~ 0.21 in 5 ks exposures. Even in a 100 ks 
exposure, the maximum redshift for temperature determi- 
nation for T — 2 keV clusters stretches only to z ~ 0.72. 
(Based on the AF98 Lx — T relation, we expect aT = 2 
keV cluster to have a luminosity of ~ 0.7 x 10 44 erg s _1 .) 
By contrast, hotter clusters, which are brighter, yield tem- 
perature measurements to higher redshifts. For example, 
temperatures will be measured for Lx > L+ (T > 5.5 keV) 
clusters to z ~ 0.6 (1.1, 2.0) in 5 (22.3, 100) ks exposures. 

5.3. Expected Catalogue Properties 



In §5.1 and § |5.2| we computed the luminosity thresh- 
old for cluster detection (Lg, a ) and temperature estimation 
(Lt), respectively, for 144,000 different combinations of 
Qq, Q\ 7 T, z, 9 and t. Combining these luminosity thresh- 
olds with the results of §|| allows us to estimate how many 
clusters will be included in our catalogue (and for how 
many of them we can estimate a temperature). We note 
that, when doing so, we assume that clusters are randomly 
locate d on the sky, that the total areal coverage is 800 
deg 2 (§4.14) and that the pointing exposure times are dis- 
tributed according to Table ||[ The results of our catalogue 
predictions are summarized in Table |I] and Figs. and ^. 

The XCS will not have a single, well defined flux limit, 
because it will be made up of po intings with a wide dis- 
persion of exposure times (§4.14), and because - in the 
XMM band at least - count rate to flux conversion fac- 



tors are a complex function of z & T (§4.6) and 9 (§4.9). 
Despite this, we have been able to estimate an effective 
flux-limit for the survey by comparing the numbers of ex- 
pected z > cluster detections (as listed on the first line 
of Table f§ 8300, 750, 61 etc.) with the N(f) values in 
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Fig. n] (after appropriate scaling from in steradians to 800 
deg 2 ). Doing so provides nine estimates of the survey flux 
limit, all of which turn out to be close to 1.5 x 10~ 14 erg 
s _1 cm -2 . Repeating the procedure for the z < 0.3 and 
z > 1.0 (by comparison with Figs. and |^ respectively) 
also yields flux limits of ~ 1.5 x 10 4 erg s _1 cm~ 2 . We 
conclude, therefore, that the effective XCS flux limit will 
be ~ 1.5 x 10~ 14 erg s _1 cm -2 , but stress that individual 
pointings in the survey will have flux limits that may be 
higher or lower than this value. 

In §^ we noted that any survey that reaches a flux limit 
of ~ 1 x 10~ 14 erg s _1 cm~ 2 will be able to detect almost 
all the T > 4 keV clusters in its survey region, irrespective 
of redshift. This is supported by Fig. ||, which shows the 
integral redshift distributions, N(z), predicted for the XCS 
in the three cosmologies we consider. The total number 
of clusters (as predicted by Press-Schechter theory) are 
depicted by solid curves, whereas the number of expected 
XCS detections are depicted by dashed curves. For T > 6 
keV clusters (top panel) the two curves are coincident out 
to z ~ 1.4, meaning that we can expect to detect all T > 6 
keV clusters at z < 1.4. From the middle panel, we can 
see that incompleteness sets in earlier for the T > 4 keV 
clusters, but even so we can expect to detect almost all 
T > 4 keV clusters out to at least z ~ 1. By contrast, we 
expect to be incomplete in terms of T > 2 keV clusters by 
z ~ 0.5 and, by z ~ 2, we can expect to be detecting only 
20% of the T > 2 keV clusters (if O = 0.3) in our survey 
region. 

Also shown on Fig. |^ are our predictions for the number 
of clusters that we will detect with sufficient signal-to- 
noise to be able to estimate temperatures (dotted curves). 
These numbers are also given in parentheses in Table [I[ 
We will obtain temperatures for all T > 6, T > 4 and 
T > 2 keV clusters out to z ~ 0.7, z ~ 0.5 and z ~ 0.3 
respectively. To further illustrate the expected proper- 
ties of the XCS, we plot in Fig. |^ the differential redshift 
distribution as a function of cosmology and temperature. 
This figure shows how many clusters will be detected in 
each Az = 0.05 bin when Oo = 1 (solid lines, square sym- 
bols), fio = 0.3, Oa = (dashed lines, star symbols) and 
O = 0.3, Oa = 0.7 (dotted lines, circular symbols). In or- 
der to differentiate between curves representing the num- 
ber of detections and curves representing the number of 
clusters with temperature estimates, we have used solid 
and open symbols respectively. 

6. DISCUSSION AND CONCLUSIONS 

A cluster catalogue of the quality described in § pT^ 
would have a great many uses. Here we describe briefly 
a subset of these, to give a feel for the kind of science 
which the XCS will make possible. We also discuss some 
caveats relating to the methods used herein, and present 
our conclusions. 

6.1. Science from the Catalogue 

The science that can be derived from the XCS can be 
loosely divided into two categories: science which can be 
obtained directly from the catalogue itself (for the most 
part assuming that follow-up observations have provided 
cluster redshifts and enabled temperature determination 
where possible), and those future projects which can build 
upon the XCS data. We give a few examples of both sorts 



of project here. 

6.1.1. Constraints on Cosmological Parameters 

The XCS cluster temperature and redshift distributions 
can be used as a direct probe of the cosmological param- 
eters O and Qa- The survey's size, redshift distribution 
and selection criteria are ideally suited to this task. The 
XCS will provide stringent constraints on Oo and has the 
potential to offer the first constraint on Oa from clus- 
ter number density evolution. The power of the XCS 
to constrain these parameters is clearly demonstrated in 
Fig. |5|, from which it is apparent that there is about an or- 
der of magnitude difference between the number of high- 
temperature (T > 4 keV) clusters in the Oq = 1 case 
compared to either of the two Oo = 0.3 cases. At lower 
temperatures (T > 2 keV) the differences between the 
various models are less apparent, demonstrating that it is 
important to concentrate on the high-temperature systems 
when attempting to measure cosmological parameters. 

It is beyond the scope of this paper to make detailed 
predictions of the errors on Oo that would result from the 
XCS. To do so would require a careful tracking of the theo- 
retical uncertainties in the number density predictions, es- 
pecially those connected with the amplitude of the power 
spectrum, and ideally would also take into account the 
weak sensitivity of the predictions to quantities such as 
the power spectrum shape. Further, the modeling of the 
observational errors, and in particular the cosmic variance 
contribution - which assesses the extent to which the ob- 
servations might be a statistical fluke - is a subtle business 
requiring detailed Monte-Carlo and probably TV-body sim- 
ulations. The former can only be carried out in detail once 
the true distribution of observing times, and the fraction 
of usable pointings, is known. It may also prove necessary 
t o mod el evolution in the temperature-luminosity relation 
(§ 16.1.21 ) . (It is for these reasons, we are unable to add error 
bars to the predictions in Figures [l], |[ [| || and 0.) 

In order to go beyond measurements of Oo and start to 
constrain Oa, one must study the z > 1 population; from 
Fig. ^ we can see that there is little difference between 
the cluster number density evolution predictions for the 
two Oo = 0.3 cosmologies below z = 1. But, for z > 1, 
the number density of galaxy clusters for Oo = 0.3 in open 
models is more than twice that in flat models. Once all the 
z > 1 clusters detected have measured redshifts and tem- 
peratures, it should be possible to constrain Oa- However, 
in view of the modeling uncertainties described above, it 
is premature to try and assess how well that can be done, 
as this will only become apparent when the actual data 
are available. 

It is important to note that the cosmological constraints 
derived from the XCS will be important even in the era 
of sensitive CMB anisotropy experiments such as Planck, 
because the cluster measurements can help to break de- 
generacies in cosmological parameters inherent in CMB 
analyses. Since the microwave anisotropy is expressed in 
terms of angular scales on the sky, the cosmological param- 
eters Oo and Oa are only constrained in the combination in 
which they arise in the angular diameter distance at the 
redshift of the last scattering surface — i.e. such obser- 
vations can only constrain the Universe to lie somewhere 
along a line in the (Oq,Oa) parameter space, rather than 
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at a single point. This degeneracy can be broken, in prin- 
ciple, in two ways. On the very largest angular scales, it is 
weakly broken by the integrated part of the Sachs- Wolfe 
effect, but these low-order multipoles suffer a large 'cosmic 
variance' which limits the accuracy with which they can 
be estimated. On small angular scales, it is mildly bro- 
ken by gravitational lensing effects. Both these effects are 
small, and it is expected that the degeneracy will largely 
still be present even after Planck has flown (Bond, Efs- 
tathiou & Tegmark 1997). One therefore requires other 
types of observations to break the degeneracy. (See Bah- 
call et al. 1999 for an overview of how different methods 
can be used in combination to constrain cosmological pa- 
rameters.) The sorts of measurements that can be used to 
break the CMB fio~^A degeneracy include cluster number 
density evolution, large-scale structure analyses, and the 
magnitude-redshift relation for Type la supernovae (Perl- 
mutter et al. 1999). Because each of these methods has 
intrinsic biases, it is important to pursue all of them, to as 
high an accuracy as possible, in order to derive a consistent 
model of the Universe. 

6.1.2. Evolution of X-ray Properties 

The XCS will greatly improve our understanding of how 
cluster properties, such as luminosity, temperature, metal- 
licity, gas mass fraction, core radius etc., evolve. Impor- 
tantly, it will allow the luminosity-temperature (Lx — T) 
relation to be measured in a coherent fashion over a wide 
redshift range. The quantification of Lx — T evolution is 
crucial to our ambition to measure cosmologicalparame- 
ters from cluster number densities. From Table m we can 
sec that the XCS will yield temperature measurements for 
1800 clusters with T > 2 keV in the low-density flat cos- 
mology. These measurements will yield the most accurate 
derivation of the Lx — T relation to date. Not only will 
the derivation include a great many more clusters, but, for 
the first time, these clusters will have been drawn from a 
single, statistically complete sample. 

6.1.3. Gravitational Lensing by XCS Clusters 

X-ray clusters magnify background galaxies via gravi- 
tational lensing. This effect is well known at optical wave- 
lengths (e.g Luppino et al. 1999), but, as shown by Smail, 
Ivison & Blain (1997), it is particularly exciting in the sub- 
millimetcr. Here, the combination of the lensing amplifi- 
cation and the positive K-correction in the submillimeter 
(resulting from the sharp decline in the spectral energy 
distribution of starburst galaxies longward of ~ 100/im) 
means such galaxies can be readily detected to extremely 
high redshift (z > 5). The follow-up of lensed galax- 
ies around XCS clusters with the coming generation of 
(sub)millimeter instruments, such as the Atacama Large 
Millimeter Array, will provide an important insight into 
the star formation history of the Universe. Lensing sig- 
nals can also be used to measure total masses and mass 
profiles of clusters in a manner complementary to X-ray 
methods (e.g. Squires et al. 1997). 

6.1.4. Sunyaev-Zel'dovich Follow-Up 

The Sunyaev-Zel'dovich (SZ) effect (Sunyaev & 
Zel'dovich 1972) describes the inverse Compton scattering 
of cosmic microwave background (CMB) photons to higher 



energies via interactions with hot electrons in the ICM. 
Measurements of the SZ effect, in combination with X-ray 
observations, provide a useful cosmological tool. They can 
be used to constrain the value of the Hubble Parameter 
(Ho, e.g. Birkinshaw 1999), the universal baryon fraction 
(e.g. Grego 1998), cluster peculiar velocities (e.g. Holzapfel 
et al. 1997) and have the potential to place powerful con- 
straints on the value of f2o (e.g. Bartlett et al. 1998). 
Therefore, SZ follow-up of XCS clusters will yield many 
important results, the most obvious of which would be the 
measurement of Hq as a function of redshift. This Hq(z) 
measurement would take advantage of the large number of 
high z clusters in the XCS and of the fact that the SZ effect 
is redshift independent. The XCS also has the potential to 
provide the required X-ray follow-up for blind SZ-surveys 
(such as that proposed by Carlstrom et al. 2000). These 
blind surveys hope to take advantage of the redshift inde- 
pendence of the SZ effect in order to detect very distant 
clusters. 

6.1.5. Analysis of CMB Foregrounds 

The limit to which the MAP and Planck satellites can 
determine the power spectrum of CMB anisotropies on 
small scales is likely to be set by the effectiveness of the 
foreground analyses. One of the major sources of CMB 
foreground confusion is the SZ signal from X-ray clusters 
of galaxies. Since the SZ effect is approximately redshift 
independent, clusters at all distances will contribute to the 
foreground signal. In a low-density cosmology, the mean 
SZ signal comes from a broad range of redshifts out to 
z ~ 2 (da Silva et al. 2000). The XCS will play a cru- 
cial role in the understanding of this signal, because it will 
provide a statistical description of the cluster population 
out to high redshifts. Moreover, in the regions covered by 
the XCS, it will be possible to mask out the signal from 
individual clusters from the CMB maps. 

6.2. Limitations of our Calculations 

As detailed in §|[ our simulations rest on a set of as- 
sumptions and simplifications, which were chosen because 
they all seem reasonable given current knowledge. There 
are, however, some limitations to the accuracy with which 
we can predict expectations for the XCS, and we discuss 
some of them below. 

6.2.1. Contamination by Low-Mass Groups 

The application of an extent criterion will not be suffi- 
cient to remove all the contamination in the cluster can- 
didate list. Low-mass groups (including "fossil groups"; 
Ponman et al. 1994; Vikhlinin et al. 1999; Romer et 
al. 2000), and some very low-redshift galaxies, will also en- 
ter the list by virtue of their extent. Low mass (and hence 
low temperature) groups are certainly interesting objects; 
they provide invaluable insight into the processes of el- 
liptical galaxy evolution, metal enrichment in the intra- 
cluster medium, and the dynamics of extended dark halos 
(Mulchaey & Zabludoff 1998). However, they have a very 
limited role to play in the derivation of Qq and f^A (be- 
cause of the increasing degeneracy between models as the 
temperature limit is decreased, see Fig. |). The Press- 
Schechter formalism becomes unreliable below T ~ 2 keV, 
so we are not able to predict the number of T < 2 keV 
groups that will be detected by the XCS. For typical values 
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of the temperature and bolometric luminosity of the intra- 
group medium (T = 1 keV, L x = 10 42 erg s _1 , Mulchaey 
& Zabludoff 1998), we have estimated the maximum red- 
shift at which a group would be detected by our survey 
to be z ~ 0.05 (0.09, 0.17) in 5 (22.3, 100) ks. We are 
actively investigating ways to flag potential T < 2 keV 
objects using a combination of extent, XMM spectra and 
cross-correlations with optical sky survey data. 

6.2.2. Contamination by Point Source Emission 



We assume in §5.3 that every cluster we detect at > 8er 



will be included in the final XCS cluster catalogue, but this 
may not always be the case. Some clusters might be ex- 
cluded if they are contaminated by point source emission, 
which might originate from an active galaxy inside the 
cluster, a foreground object, such as an M star, or a back- 
ground object, such as a quasar. Romer et al. (2000) de- 
scribe the case of an extended X-ray source (RXJ0947.8) 
which was excluded from their cluster catalogue, on the 
grounds of it being coincident with a z — 0.63 quasar, 
despite there being a spectroscopically confirmed cluster 
at the same position and redshift. A total of four clus- 
ters were rejected from the Bright SHARC cluster cat- 
alogue because the quality of the ROSAT data did not 
permit the cluster flux to be disentangled from that of a 
contaminating point source. Romer et al. (2000) claim 
that two of these systems probably have sufficient, uncon- 
taminated, flux to merit inclusion in the Bright SHARC 
cluster catalogue, which corresponds to an incompleteness 
of the whole catalogue at the 5% level. We expect the 
incompleteness level to be much lower than this for the 
XCS since the improved spatial resolution of XMM over 
ROSAT will significantly enhance our ability to mask out 
point source contamination when measuring cluster fluxes. 
We expect therefore that the wrongful exclusion of clusters 
from the XCS catalogue will occur only very rarely and 
have an insignificant effect on our ability to use the XCS 
as a cosmological tool. 

6.2.3. Effect of Assumptions about the Cluster Model 

Perhaps the greatest uncertainty in our calculations 
comes from the simplified model of the distribution and 
state of the intracluster medium we employed. Our use of 
the isothermal /3-model was justified in part by the results 
of Mohr et al. (1999), who showed that it well described the 
azimuthally-averaged properties of known clusters. How- 
ever, this work was carried out in the ROSAT bandpass 
(0.5-2.0 keV), over which the emissivity of the X-ray gas 
is almost insensitive to cluster temperature for T ;> 2 keV. 
And, as emphasized recently by Ettori (1999), the assump- 
tion of a simple isothermal /3-model will lead to significant 
errors when a cluster with a significant temperature gra- 
dient is observed in a broad band bracketing the energy 
corresponding to its mean temperature. Evidence to sug- 
gest that the cluster gas is not isothermal comes from spa- 
tially resolved cluster temperature maps (e.g. Markevitch 
et al. 1998) and from the so-called "/3 discrepancy" (e.g. 
Sarazin 1988; Bahcall & Lubin 1993), which describes the 
fact that fitted values of /3 are not consistent with the val- 
ues expected from the combination of cluster temperatures 
with galaxy velocity dispersions. 

We have also ignored the effect of cooling flows in the 
cluster core. A significant fraction of relaxed clusters have 



regions of cool, dense gas in their cores (e.g. AF98) and, 
as Ettori (2000) has pointed out, a modified version of 
equation would be more appropriate to describe such 
clusters. 

More fundamentally, it is possible, of course, that the 
clusters we detect at high redshift will not be virialised 
systems. Clusters in the process of formation may have 
significant non-thermal components to their X-ray lumi- 
nosities, for example from shocks resulting from subcluster 
merging. A classic example of such a cluster is RX J0152.7 
(Romer et al. 2000; Ebeling et al. 2000), which has a high 
total luminosity (8.26 x 10 44 erg s _1 , z = 0.83) but is made 
up of at least 2 components. It is not possible to predict, 
at this stage, what the net effect of unvirialised systems 
will be on the properties of the XCS, but the spatial and 
spectral resolution of EPIC should help us to recognize 
such systems. The XCS may even show that clusters are 
not suitable as cosmological probes above a certain red- 
shift: indeed, perhaps effects such as these lie behind the 
detection to date of a (possibly) surprisingly large num- 
ber of massive clusters at high redshift (Luppino & Gioia 
1995; Donahue 1996; Luppino & Kaiser 1997; Donahue et 
al. 1998; Eke et al. 1998), which has been claimed to be 
troublesome for conventional models of structure forma- 
tion. 

6.3. Conclusions 

We have predicted the expected properties of a serendip- 
itous cluster survey based on archival XMM pointing data. 
We have done this using simulations which combine a theo- 
retical model of the properties of the cluster population, as 
a function of cosmology, with a detailed description of the 
characteristics of the EPIC camera, and a generic model 
for cluster surface brightness profiles. We have shown that 
the catalogue that would result from such a survey will 
surpass existing catalogues of high-redshift (z > 0.3) clus- 
ters in both size, quality and redshift coverage, while, at 
low redshifts (z < 0.3) the catalogue will yield many more 
cluster temperature measurements than have ever been 
measured before. 

It is clear that, while the methods presented here may be 
adequate to yield reasonably realistic predictions for what 
we can expect to get from the XCS, the actual analysis 
of the data from the survey will require a more sophis- 
ticated approach, informed by detailed physical models 
resulting from pointed observations of individual clusters 
made by Chandra and XMM itself. This must, however, 
be tempered by the requirement that the final set of clus- 
ter selection criteria be readily modeled by Monte-Carlo 
methods: the XCS will have its greatest impact in sta- 
tistical analyses, so it must be constructed in such a way 
that its selection function can be well understood. 
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Fig. 1. — The expected number of galaxy clusters across the whole sky with X-ray temperatures in excess of 6 keV (upper panel), 4 keV 
(middle panel) and 2 keV (lower panel), as a function of the catalogue flux threshold in the 0.5-2.0 keV band. The Qo = 1 case is the solid 
line, while for Qq = 0.3 the flat case is shown as dashed and the open case as dotted. 
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Fig. 2. — As Fig. [ij but showing only clusters with z > 1. 
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Fig. 3. — As Fig. Ill but showing only clusters with z < 0.3. 
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Fig. 4. — Fitted temperatures versus input temperatures for 4 different combinations of redshift and background contamination. All 
spectra were created using f akeit in XSPEC. The dotted lines show input temperature plus (upper line) and minus (lower line) 20%. 
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Fig. 5. — The cumulative redshift distribution, JV(< z), of galaxy clusters per 800 square degrees with X-ray temperature in excess of 6 
KcV (upper panel), 4 keV (middle panel) and 2 keV (lower panel). The solid lines show the result one would obtain if there was no limitation 
on the detectable flux. The dashed and dotted lines show our predictions for the XCS, where the dashed line represents the expected number 
of > 8cr detections and the dotted line represents the expected number of clusters bright enough to allow temperature measurements. 
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Fig. 6. — The predicted redshift distributions (evaluated in bins of width Az = 0.05) corresponding to the cumulative counts of Fig. |5j. The 
filled and empty symbols denote clusters detected and with temperatures estimated, respectively, in the three cosmological models: Fjnstein 
- de Sitter (solid line), low density open (dashed) and low density flat (dotted). 



Table 1 

Coefficients of the quadratic fits to K-corrections for Raymond-Smith spectra. 
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0.5-10.0 keV 
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Table 2 

Conversion factors between observed and pseudo-bolometric luminosity for Raymond-Smith spectra. 



Temperature (keV) 


0.5-10.0 keV 


0.5-2.0 keV 


1 


1.75287 


2.00900 


2 


1.48716 


2.37600 


3 


1.36920 


2.69564 


4 


1.33470 


3.04131 


5 


1.33501 


3.38720 


6 


1.35271 


3.72926 
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1.38674 


4.04849 
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4.36129 
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4.66436 
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4.95957 


11 
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Table 3 

Distribution of XMM exposure times in the GTO and AOl observing cycles. 



time (ks) 540 1045 i!T20 2CT25 2541) 30455 35^40 4CT45 45^50 5CT55 55-100 

percentage (GTO) 23.16 18.03 7.63 17.50 4.87 6.58 1.45 7.24 1.58 6.58 5.39 
percentage (API) 23.79 18.97 4.41 18.81 5.56 9.08 1.31 4.66 0.57 7.23 5.56 
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Table 4 

The expected number of clusters detected in an XMM serendipitous survey covering 800 deg 2 , for three 

DIFFERENT COSMOLOGICAL MODELS. THE MAIN NUMBERS ARE FOR DETECTIONS, WHILE THE NUMBERS IN PARENTHESES ARE 

DETECTIONS WITH SUFFICIENT FLUX TO YIELD TEMPERATURES. 





Q = 0.3, A = 0.7 

T>2 T>4 T>6 


z > 
z > 0.3 
z > 1 


8300 (1800) 750 (320) 61 (42) 
7600 (1200) 700 (270) 54 (36) 
750 (6) 170 (6) 12 (2) 




fi = 1.0, fU = 0.0 
T>2 T>4 T>6 


z > 
z > 0.3 
z > 1 


2600 (1200) 80 (70) 5 (5) 
1900 (570) 50 (40) 2 (2) 
46 (0) 1 (0) (0) 




O = 0.3, O a = 0.0 
T>2 T>4 T>6 


z > 
z > 0.3 
z > 1 


9700 (2000) 1100 (380) 110 (56) 
9000 (1400) 1100 (330) 100 (51) 
1700 (26) 480 (24) 50 (9) 



